function rest = L2error(uex,u)

rest=sqrt(sum((uex-u).^2)*1/(size(uex,1)+1));


